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ABSTRACT 


This  report  describes  the  determination  of  empirical  parameters  for  a  fracture  model  which  can  be 
used  to  predict  the  uniaxial  tensile  strength  of  composite  laminates  containing  notches  of  various 
sizes.  The  parameters  required  are  the  characteristic  dimension  and  the  unnotched  tensile 
strength  of  the  subject  laminate.  A  data  base  including  notched  and  unnotched  strengths  of  the 
subject  laminates  and  elastic  ply  properties  was  generated  for  the  determination  of  these 
parameters  using  a  statistical  procedure  based  on  the  two  parameter  Weibull  distribution.  Three 
second  generation  high  toughness  graphite  fibre  reinforced  composites,  IM6/5245C,  IM6/F584 
and  IM6/1806,  and  two  first  generation  brittle  graphite/epoxy  composites,  AS4/3501-6  and 
T300/5208,  were  used  to  fabricate  the  specimens  containing  either  a  through-the-thickness 
cylindrical  hole  or  a  countersunk  hole. 


The  two  parameter  fracture  model  assessed  was  found  to  be  capable  of  predicting  accurately  the 
tensile  notched  strength  of  composite  laminates.  The  characteristic  dimension  was  found  to  be 
dependent  on  a  variety  of  intrinsic  and  extrinsic  factors;  therefore,  it  must  be  determined  by 
experiment  for  each  new  composite  material,  laminate  configuration,  etc.  Results  showed  that  a 
decrease  in  the  characteristic  dimension  indicated  an  increase  in  notch  sensitivity.  Since  the 
characteristic  dimensions  of  the  three  second  generation  composites  were  found  to  be  larger 
than  those  of  the  first  generation  composites,  they  were  less  notch  sensitive.  The  notched 
strength  of  a  countersunk  hole  was  predicted  by  approximating  the  countersunk  hole  by  a  uniform 
hole  with  a  diameter  equal  to  the  outer  diameter  of  the  countersunk  hole. 


RESUME 

Le  present  rapport  decrit  la  determination  de  parametres  empiriques  en  vue  de  I’elaboration  d'un 
modele  de  rupture  permettant  de  predire  la  resistance  a  la  traction  uniaxiale  de  materiaux 
composites  stratifies  comportant  des  entailles  de  differentes  failles.  Les  parametres  requis  sont  la 
dimension  caracteristique  et  la  resistance  a  la  traction  en  I'absence  d’entailles  du  materiau  stratifie. 
Une  base  de  donnees  contenant  les  resistances  des  materiaux  stratifies  concernes  avec  et  sans 
entailles  ainsi  que  les  proprietes  elastiques  des  plis  a  ete  etablie;  elle  permet  de  determiner  ces 
parametres  a  I'aide  d'une  methode  statistique  basee  sur  la  distribution  de  Weibull  des  deux 
parametres.  Trois  materiaux  composites  renforces  de  fibres  de  graphite  tres  resistants  aux  chocs, 
de  deuxieme  generation,  IM6/5245C,  1M6/F584  et  IM6/1806,  et  deux  materiaux  composites 
graphite/resine  epoxy  cassants  de  premieres  generation,  AS4/3501-6  et  T300/5208,  ont  ete 
utilises  pour  la  fabrication  d'eprouvettes  contenant  soit  un  trou  cylindrique  traversant  toute 
I’epaisseur,  soit  une  fraisure. 

Le  modele  de  rupture  a  deux  parametres  evalue  s'est  avere  capable  de  predire  avec  precision  la 
resistance  a  la  traction  de  materiaux  composites  stratifies  comportant  des  entailles.  On  a  constate 
que  la  dimension  caracteristique  dependait  de  divers  facteurs  intrinseques  et  extrinseques:  elle 
doit  done  etre  determinee  experimentalement  pour  cheque  nouveau  materiau  composite, 
cheque  type  de  stratification,  etc.  Les  resultats  monirent  qu'une  reduction  de  la  dimension 
caracteristique  s'accompagne  d'une  augmentation  de  la  sensibilite  a  I'effet  d'entaille.  Etant  donne 
que  les  dimensions  caracteristiques  des  trois  materiaux  composites  de  deuxieme  generation 
etaient  superieures  a  celles  des  trois  materiaux  composites  de  premiere  generation,  les  materiaux 
de  deuxieme  generation  etaient  moins  sensibles  a  I'effet  d'entaille.  Le  resistance  a  la  traction  des 
materiaux  comportant  une  fraisure  etait  calculee  en  assimilant  la  fraisure  a  un  trou  uniforme  de 
diametre  egal  au  di=‘motr?  exterieur  de  la  fraisure. 

i  i  i 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  .  iii 

LIST  OF  TABLES .  v 

LIST  OF  FIGURES .  vi 

LIST  OF  ABBREVIATIONS .  vii 

1.0  INTRODUCTION .  1 

2.0  NOTCHED  STRENGTH  MODEL .  2 

3.0  EXPERIMENTAL  PROCEDURES  AND  RESULTS  .  4 

3.1  Unnotched  Tensile  Tests .  4 

3.2  Notched  Tensile  Tests .  5 

3.3  Data  Reduction .  5 

3.4  Determination  of  Parameters  for  the  Fracture  Model .  6 

3.4.1  Statistical  analysis  of  unnotched  strengths  .  7 

3.4.2  Statistical  analysis  of  notched  strengths  .  7 

3.4.3  Determination  of  a^  .  8 

4.0  ASSESSMENT  OF  FRACTURE  MODEL  .  8 

4.1  Notched  Strength  Prediction  .  9 

5.0  CONCLUSIONS .  11 

6.0  ACKNOWLEDGEMENT .  12 

7.0  REFERENCES .  12 


IV 


LIST  OF  TABLES 


1.  Mechanical  Properties  of  Composite  Materials 

2.  Stress  Concentration  Factors  for  Composite  Laminates 

3.  Unnotched  Tensile  Test  Results 

4.  Test  Matrix  for  Notched  Specimens 

5.  Notched  Tensile  Test  Results 

6.  Results  of  the  Determination  of  the  Weibull  Parameters  for  Unnotched  Test  Data 

7.  Results  of  the  Determination  of  the  Weibull  Parameters  for  Notched  Test  Data 

8.  Characteristic  Dimensions  for  IM6/5245C  and  AS4/3501-6 


V 


LIST  OF  FIGURES 


1.  Average  Stress  Failure  Criterion  for  a  Circular  Through-the-Thickness  Hole 

2.  Effect  of  the  Characteristic  Dimension,  a^,,  on  the  Notched  Strength  Prediction  for 
Laminates  Containing  a  Circular  Hole  Having:  (a)  Kj°°  =  3.0,  and  (b)  K]“  =  6.0 

3.  A  Typical  Stress  vs  Strain  Plot  for  an  Unnotched  Tensile  Test 

4.  A  Typical  Setup  for  a  Notched  Tensile  Test 

5.  A  Typical  Stress  vs  Strain  Plot  for  a  Notched  Tensile  Test 

6.  Weibull  Probability  Plot  for  IM6/5245C  Unnotched  Strength  Data 

7.  Weibull  Probability  Plot  for  AS4/3501-6  Unnotched  Strength  Data 

8.  Weibull  Probability  Plot  for  IM6/5245C  Notched  Strength  Data 

9.  Weibull  Probability  Plot  for  AS4/3501-6  Notched  Strength  Data 

10.  Comparison  between  Predicted  /a ^  and  Test  Results  for  IM6/5245C 

11.  Comparison  between  ja ^  and  Test  Results  for  Different  Composites 

12.  Effects  of  Characteristic  Dimension  uf  Notch  Sensitivity 

13.  Effects  on  Failure  Mode  Resulting  from  Increasing  Hole  Size 

14.  Flat  Fracture  of  an  IM6/5245C  Laminate  Containing  a  Countersunk  Hole 


VI 


LIST  OF  ABBREVIATIONS 


So 

A, 

D 


E 


ri’ 


E  ,  E 

X’  y 

Gi2 


K, 


L 


R 

W 


X., 

Xt2 

Xc, 

Xc:2 


X. 


B 


■-x* 

N 


characteristic  dimension 
in-plane  laminate  stiffness 
hole  diameter 

moduli  of  elasticity  of  lamina  in  0°  and  90“  directions 
effective  moduli  of  elasticity  of  laminate  in  x  and  y  directions 
shear  modulus  of  lamina 
effective  shear  modulus  of  laminate 

orthotropic  stress  concentration  factor  for  finite  width  laminate 

orthotropic  stress  concentration  factor  for  infinite  width  laminate 

rank  number 

laminate  length 

hole  radius 

laminate  width 

average  tensile  strength  of  lamina  in  0“  direction 
average  tensile  stiength  of  lamina  in  90°  direction 
average  compressive  strength  of  iamina  in  0°  direction 
average  compressive  strength  of  lamina  in  90°  direction 
abscissa  of  the  Weibull  probability  plot 
ordinate  of  the  Weibull  probability  plot 

Weibull  unnotched  characteristic  strength  parameter  of  laminate 
Weibull  notched  characteristic  strength  parameter  o*'  infinite  width  laminate 

vii 


1.0  INTRODUCTION 


In  recent  years,  extensive  research  and  development  has  been  conducted  with  the  aim 
of  increasing  the  application  of  composite  materials  in  aircraft,  particularly  in  the  primary 
airframe  structures.  A  major  area  of  study  involves  the  notch  sensitivity  of  composite 
laminates.  The  results  of  these  investigations  have  produced  a  variety  of  analytical 
techniques  ranging  from  comprehensive  numerical  methods  to  simplified  fracture  models 
that  are  capable  of  addressing  the  complexity  of  failure  processes  and  damage  progression 
in  notched  composite  laminates. 

Several  fracture  models  based  on  linear  elastic  fracture  mechanics  (LEF.M)  have 
been  proposed  extending  their  applications  from  metals  to  composites  containing  notches 
Wu  indicated  that  LEFM  can  only  be  applied  directly  to  composites  when  very 
specific  conditions  are  satisfied.  The  application  of  LEFM  has  been  restricted  because  in 
the  majority  of  cases  no  self-similar  damage  growth  is  observed.  The  complexity  of  the 
damage  mechanism  at  the  notch  also  raises  serious  questions  as  to  the  applicability  of 
LEFM  to  composites.  On  the  microscopic  level  the  failure  modes  occurring  within  the 
damage  zone  at  the  notch  appear  in  the  form  of  fiber  pull-out,  matrix  micro-cracking  and 
fiber-matrix  interfacial  failure.  On  the  macroscopic  level  the  major  failure  modes  are 
delamination,  matrix  cracking  along  the  fiber  direction  and  fiber  breakage. 

It  is  very  difficult  to  determine  accurately  the  stress  distribution  and  to  account  for 
the  variety  of  failure  modes  at  the  notch.  As  a  result  semi-empirical  fracture  models 
which  ignore  the  physical  state  of  the  damage  zone,  have  been  developed.  These  fracture 
models  are  based  on  elastic  stress  or  strain  distribution  at  a  distance  from  the  notch  and 
require  the  use  of  experimentally  determined  parameters  for  the  prediction  of  the  notched 
strength  of  composites.  The  fracture  behaviour  of  notched  composite  laminates  has  been 
shown  to  be  affected  by  many  intrinsic  and  extrinsic  variables.  Intrinsic  variables  such  as 
laminate  configuration,  stacking  sequence,  constituent  properties,  fiber  volume  fraction, 
fiber-matrix  interface  characteristics  and  fabrication  procedure,  and  extrinsic  variables  such 
as  specimen  geometry,  notch  configuration,  moisture  content,  temperature  and  loading  rate 
are  known  to  have  significant  influence  on  the  notched  strength  and  the  failure  mode  of 
composite  laminates  Therefore,  the  effect  of  these  variables  on  the  parameters  of  the 
fracture  model  must  be  empirically  determined. 

Numerous  experimental  investigations  have  been  conducted  in  the  past  fifteen  years 
with  the  aim  of  deriving  the  parameters  required  by  the  fracture  models  to  account  for  the 
effect  of  the  varous  intrinsic  and  extrinsic  variables.  The  majority  of  the  investigations  dealt 
with  the  notched  strength  of  first  generation  graphite,''epoxy  composites  subjected  to  quasi¬ 
static  uniaxial  tensile  loading  in  ambient  temperature  and  moisture  conditions.  Ver>  limited 
data  are  available  in  the  literature  for  boron/aluminium,  graphitc/polyimide  and  second 
generation  high  toughness  composites. 


The  objectives  of  lliis  investigation  were  to  establish  a  data  base  for  assessing  the 
tensile  notch  sensitivity  of  several  second  generation  high  toughness  composites  and  to 
establish  the  parameters  of  a  fracture  model  for  the  prediction  of  notched  strength  of 
composite  laminates  containing  a  through-the-thickncss  hole.  Both  countersunk  and 
noncountersunk  holes  were  investigated.  'I'he  notch  sensiti\ity  due  to  an  open  hole  and 
impact  damage  under  compressive  loading  are  reported  elsewhere  A  w  idely  used  two- 

parameter  fracture  model  developed  by  Whitney  and  Nuisn^er  was  adopted  in  this  work. 
For  the  sake  of  completeness  the  theory  of  this  model  is  outlined  in  Section  2.  Experimental 
procedures  including  testing  techniciues  and  statistical  reduction  of  data  to  derive  the 
parameters  for  notched  strength  prediction  are  presented  in  Section  3.  The  assessment  of 
the  fracture  model  by  using  it  to  perform  notched  strength  predictions  is  gi\en  in  Section 
4.  Conclusions  are  provided  in  Section  5. 


2.0  NOTCHED  STRENGTH  MODEE 


The  fracture  model  developed  by  Whitney  and  Nuismer'’'  was  used  to  predict  the 
tensile  notched  strength  of  composite  laminates  containing  a  through-the-thickness  hole.  The 
model  is  based  on  the  average  stress  criterion  which  hypothesizes  that  failure  will  occur 
when  the  stress  at  a  characteristic  dimension,  a^,.  from  the  edge  of  the  hole  reaches  the 
unnotched  strength,  a„,  of  the  laminate  (Fig.  1).  This  criterion  can  be  expressed  by: 

~  (1) 


The  model  requires  the  knowledge  of  a,,,  a,,,  and  the  normal  stress,  a^,  at  the  hole 
in  order  to  predict  the  notched  strength.  It  is  assumed  that  a„  is  a  material  parameter  and 
as  a  result  it  is  independent  of  the  size  of  the  hole.  Therefore,  this  fracture  model  requires 
only  tw  o  parameters,  a„  and  a^,,  for  the  prediction  of  tensile  strength  of  composite  laminates 
containing  notches  of  various  sizes. 

Thea„  of  the  laminate  can  be  determined  bv  performing  a  tensile  test  on  the  subject 
laminate  or  by  analysis  using  lamination  plate  theory  and  material  allowables.  The 
characteristic  dimension  must  be  determined  bv  corroborating  with  experimental  data.  In 
this  investigation,  both  the  imi'.otched  strength  and  the  charactei istic  dimension  were 
determined  experimentally,  d'he  experimental  procedures  used  will  be  described  in  Section 
3. 


l  or  an  infinite  vJi  thotropic  plate  containing  a  circular  hole  of  radius  R  and  subjected 
to  an  uniform  stress,  o  '  ,  applied  parallel  to  the  v-axis  at  infinit\  (1  ig.  1).  the  normal 
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stress,  Oy  along  the  x-axis  in  the  remaining  ligament  of  the  laminate  can  be  expressed 
approximately  by 
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where  Kj.”^  is  the  orthotropic  stress  concentration  factor  for  an  infinite  width  plate  as 
expressed  by 
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where  A-  are  the  in-plane  laminate  stiffnesses'"^’'  as  determined  from  laminated  plate  theor\ . 

The  exact  solution  for  the  in-plane  stress  distribution  in  an  infinite  orthctropic  plate 
containing  a  circular  hole  and  subjected  to  remote  uniaxial  tension  has  been  formulated  by 
Lekhnitskii  However,  the  exact  solution,  using  a  complex  variable  mapping  approach, 
is  rather  complicated.  In  order  to  facilitate  computation,  a  polynomial  approximation  of  the 
exact  solution  is  required.  Konish  and  Whitney  demonstrated  that  the  polynomial 
expression  given  in  Equation  2  provides  a  very  good  approximation  to  the  exact  solution  for 
[0],  [±45]  and  a  series  of  [0/±45]^  T300/5208  graphite/epoxy  laminates.  For  isotropic  plates 
or  laminates  with  quasi-isotropic  lay-up  (K|^’  =  3.0),  the  sixth  and  eighth  auc.  terms  in 
Equation  2  can  be  omitted. 


Combining  Equations  1  and  2  and  carrying  out  the  integration  result  in  the  ratio  of 
notched  strength  to  unnotched  strength  as  expressed  by: 
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The  lower  and  upper  limits  ot  Equation  4  are: 


(1)  For  very  large  holes:  i 


(2)  For  very  small  holes:  i 


0, 
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A  series  of  notch  sensitivity  curves  can  be  obtained  by  plotting  the  values  ofa^°°/a„ 
as  a  function  of  R.  The  notch  sensitivity  curves  demonstrating  the  effect  of  a^  on  the 
notched  strength  of  composite  laminates  containing  a  through-the-thickness  circular  hole 
with  K-p°°  =  3.0  and  Kp°°  =  6.0  are  shown  in  Fig.  2a  and  Fig.  2b  respectively. 

3.0  EXPERIMENTAL  PROCEDURES  AND  RESULTS 


The  notch  sensitivity  curves  shown  in  Fig.  2  can  be  used  to  predict  the  tensile 
notched  strength  of  a  wide  class  of  orthotropic  composite  laminates  when  and  a^,  of  the 
material  are  known.  In  this  section,  the  experimental  procedures  and  the  data  reduction 
techniques  used  to  determine  these  two  parameters  are  discussed. 

The  composite  materials  used  and  their  unidirectional  mechanical  properties  are 
given  in  Table  1.  The  unidirectional  mechanical  properties  were  obtained  from  a  data  base 
which  has  been  developed  from  the  characterization  work  on  advanced  composite  materials 
conducted  at  the  Institute  for  Aerospace  Research  (formerly  known  as  the  National 
Aeronautical  Establishment)  These  properties  were  used  to  calculate  the  in-plane 

laminate  stiffness  terms  which  were  then  substituted  into  Equation  3  to  calculate  the 
orthotropic  stress  concentration  factors  for  the  laminates  tested  in  this  work.  The  results  are 
shown  in  Table  2. 

Fabrication  of  the  composite  panels,  determination  of  the  fiber  content  using  the 
burn-off  process  and  inspection  of  the  panels  using  ultrasonic  C-scan  are  decribed  in 
References  31,  32  and  33  respectively. 

3.1  Unnotched  'Tensile  Tests 

Tensile  tests  were  performed  on  unnotched  specimens  fabricated  with  the  Narmco 
IM6/5245C,  a  second  generation  high  tougiiness  composite,  and  from  the  Hercules 
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AS4/3501-6,  a  first  generation  graphite/epoxy  composite.  The  specimens  fabricated  with  the 
IM6/5245C  were  nominally  2.5  cm  wide  by  30  cm  long  with  a  gauge  section  of  15  cm.  The 
specimens  fabricated  with  the  AS4/3501-6  were  nominally  2.5  cm  wide  by  43  cm  long  with 
a  gauge  section  of  27  cm.  These  specimen  configurations  meet  the  width  and  gai^re  length 
requirements  of  the  ASTM  Standards' About  half  of  the  specimens  tested  were  strain 
gauged.  A  pair  of  hydraulic  grips  was  used  to  hold  the  specimens  in  an  MTS  machine  with 
a  500  kN  capacity.  A  gripping  technique'^^'  involving  the  use  of  sandblasted  aluminium  shims 
was  employed  to  prevent  specimen  slippage  during  load  transfer.  The  specimen  was  loaded 
in  tension  to  failure  at  a  cross  head  displacement  rate  of  0.25  mm/min.  A  typical  stress  vs 
strain  curve  for  the  unnotched  tension  test  is  shown  in  Fig.  3.  The  results  of  the  unnotched 
tension  tests  are  summarized  in  Table  3. 

3.2  Notched  Tensile  Tests 

Tensile  tests  were  performed  on  specimens  containing  a  through-the-thickness  hole. 
Both  cylindrical  and  countersunk  holes  were  investigated.  The  specimens  tested  were 
manufactured  from  three  second  generation  high  tounghness  composites,  which  are  the 
Narmco  IM6/5245C,  the  Hexcel  IM6/F584  and  the  Cyanimide  IM6/1806,  and  from  two  first 
generation  graphite/epoxy,  which  are  the  Narmco  T300/5208  and  the  Hercules  AS4/3501-6. 
A  test  matrix  is  provided  in  Table  4. 

Most  of  the  specimens  tested  were  nominally  5  cm  wide  by  30  cm  long  with  a  gauge 
section  of  15  cm  between  hydraulic  grips.  The  specimens  fabricated  with  the  AS4/3501-6 
materials  were  nominally  3.8  cm  wide  by  45  cm  long.  The  specimens  containing  a 
countersunk  hole  were  5  cm  wide  by  42.5  cm  long.  These  specimen  configurations  resulted 
in  length-to-width  (L/W)  and  diameter-to-width  (D/W)  ratios  that  satisfy  the  requirements 
for  finite  width  correction  explained  in  Sub-section  3.3.  All  specimens  tested  were  strain 
gauged.  A  pair  of  hydraulic  grips  with  sandblasted  aluminium  shims  was  used  to  hold  the 
specimens  in  an  MTS  machine  with  a  capacity  of  500  kN.  A  typical  setup  of  the  notched 
tensile  test  is  shown  in  Fig.  4.  The  specimens  were  loaded  in  tension  to  failure  at  a  cross 
head  displacement  rate  of  0.25  mm/min.  A  typical  stress  vs  strain  curve  for  a  notched 
tension  lest  is  show  n  in  Fig.  5.  Table  5  summarizes  the  results  of  the  notched  tension  tests. 

3.3  Data  Reduction 

The  stress  at  failure  of  the  unnotched  specimens,  a„,  was  computed  based  on  the 
failure  load  and  the  average  cross  sectional  area  of  the  specimens.  The  modulus  of  the 
laminate  was  computed  by  dividing  the  stress  with  the  strain  when  the  strain  was  .004 
mm/mm.  The  stress  at  failure  of  the  notched  specimen,  Oj^-,  was  computed  by  dividing  the 
failure  load  with  the  average  gross  cross  sectional  area  of  the  specimens. 

Since  the  computed  is  for  specimens  with  a  finite  width,  a  correction  factor  is 
required  to  obtain  the  notched  strength  of  a  laminate  with  an  infinite  width,  ,  for  use 
in  the  fracture  model  (equation  4).  When  the  material  tested  is  orthotropic,  an  appropriate 
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orthotropic  finite  width  correction  factor  should  be  used.  1'he  correction  factor  can  be 
determined  accurately  using  analytical  methods  involving  finite  element  boundary 
collocation  or  boundary  integral  equation  It  should  be  noted  that  no  closed  form 
solutions  are  available  and  numerical  solutions  are  required  to  obtain  the  orthotropic 
correction  factor.  In  order  to  simplify  the  calculation,  an  isotropic  finite  width  correction 
factor''’^',  shown  in  Equation  6,  was  used  in  this  work. 


The  validity  of  employing  the  isotropic  width  correction  factor  with  orthotropic 
materials  has  been  addressed  extensively  in  the  literature.  The  use  of  the  isotropic  factor 
instead  of  the  orthotropic  factor  has  been  shown  to  be  valid  when  the  D/VV  ratio  is  less  than 
or  equal  to  1/4*^'*'  and  the  L/W  ratio  is  larger  than  3'^'*.  As  mentioned  earlier  all  specimens 
used  in  this  work  meet  these  two  requirements.  Using  Equation  6,  finite  width  corrections 
were  performed  on  the  experimentally  determined  notched  strengths  and  the  results  are 
shown  in  Table  5. 

3.4  Determination  of  Parameters  for  the  Fracture  Model 

The  fracture  model  used  in  this  investigation  has  two  unknown  parameters,  a ^  and 
a^,.  which  must  be  determined  before  it  can  be  used  to  predict  the  strength  of  notched 
laminates.  A  statistical  procedure  based  on  the  two-parameter  Weibull  distribution  was  used 
to  determine  these  two  parameters  for  the  Narmco  1M6/5245C  and  the  Hercules  AS4/3501- 
6.  The  probability  density  function  for  the  Weibull  distribution  is  expressed  by: 


/:  (  o, 


where  Of  =  stress  at  failure 

B  =  characteristic  strength  parameter 
k  =  slope  parameter 
(p  =  location  parameter 


When  0=0.  Equation  7  reduces  to  the  two-parameter  Weibull  distribution,  riie 
characteristic  strength  parameter  occurs  at  the  63.2%  failure  point  for  the  sample 
population.  The  Weibull  distribution  is  a  general  function  which  encompasses  several  well 
known  distributions.  For  example,  the  function  given  in  Equation  7  is  a  simple  exponential 
distribution  function  when  k  =  1;  the  Rayleigh  distribution  function  when  k  =  2;  and  a  good 
approximation  of  the  Normal  distribution  function  when  k  =  3..S7.  that  is,  when  the  mean 


and  the  median  values  are  equal. 

3.4.1  Statistical  analysis  of  unnotched  strengths 


The  cumulative  form  of  the  probability  density  function  shown  in  Equation  7  was 
used  to  determine  the  two  Weibull  parameters,  and  for  the  unnotched  strength  data 
shown  in  Table  3.  This  cumulative  function  is  expressed  by: 


P  (  )  -  1  -  exp 

V 


(8) 


The  unnotched  strength  data  were  first  ranked  from  minimium  to  maximum  strength. 
The  stress  at  failure  of  each  unnotched  specimen  was  then  assigned  a  rank  number,  i,  from 

I  through  n,  where  n  is  the  sample  size.  In  this  report,  P(Op)  is  estimated  by’^^': 

(9) 

n 

The  unknown  parameters,  and  k^,.  were  determined  using  a  plotting  technique. 
In  the  Weibull  probability  plot,  the  abscissa  and  the  ordinate  are  shown  in  Equations  10  and 

II  respectively. 


In 


In 


1  - 


(10) 


-  In  [  a  J 


(11) 


The  Weibull  plots  of  the  unnotched  strength  data  shown  in  Table  3  are  shown  in 
Figs.  6  and  7.  Using  linear  regression,  a  straight  line  was  fitted  through  the  data.  The 
parameter  k^,  was  obtained  by  determining  the  slope  of  the  straight  line.  The  parameter  B^, 
was  determined  by  reading  off  the  strength  value  corresponding  to  the  intersection  of  the 
fitted  line  and  a  horizontal  line  corresponding  to  P(o„)  =  63.2%.  The  results  of  the 
calculation  are  summarized  in  Table  6. 

3.4.2  Statistical  analysis  of  notched  strengths 

The  notched  strength  data  shown  in  Table  5  were  analysed  using  the  same  procedure 
discussed  in  Sub-section  3.4.1  for  the  Narmco  IM6/5245C  and  the  Hercules  AS4/3501-6 
materials.  The  Weibull  plots  of  the  notched  strength  data  arc  shown  in  Figs.  8  and  9.  'fhe 
parameters,  Bv^”^  and  kv,_,  for  the  two  materials  arc  shown  in  fable  7. 
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3.4.3  Determination  of 

The  characteristic  strength  parameters,  and  were  substituted  into  Equation 
4  to  yield  Equation  12  which  was  solved  iteratively  for  a^,.  The  results  are  shown  in  Table 
8. 


^  ^  _ 2  (  1-n _ 

Po  2  -  e  -  (  K;  -  3)  i  -  e  ) 


(12) 


4.0  ASSESSMENT  OE  FRACTURE  MODEL 


As  mentioned  earlier,  the  notch  sensitivity  of  composite  laminates  is  affected  by  many 
intrinsic  and  extrinsic  factors.  The  applicability  of  the  fracture  model  would  be  greatly 
enhanced  if  the  parameters  of  the  model  could  be  applied  universally  to  predict  the  notched 
strength  of  laminates  manufactured  with  different  composite  materials  and  laminate 
configurations  and  containing  different  notch  shapes.  A  major  requirement  for  the  fracture 
model  to  be  universal  is  that  the  characteristic  dimension  is  a  material  constant  and  is 
independent  of  specimen  geometry  and  laminate  configuration. 

Whitney  and  Nuismer  claimed  and  substantiated  with  experimental  evidence  that 
using  the  average  stress  failure  criterion  with  a  constant  value  of  a^  =  2.3  mm  resulted  in 
acceptable  predictions  for  the  fracture  strength  of  a  number  of  different  laminates  of  an 
AS/3501-5  graphite  epoxy  material  system.  Their  results  also  showed  that  using  the  model 
resulted  in  good  strength  predictions  for  countersunk  and  non-countersunk  holes,  and  for 
normal  and  slant  cracks. 

Several  investigators  attempted  to  use  the  same  fracture  model  that  Whitney  and 
Nuismer  used  for  predicting  the  notched  strength  of  different  laminates  and  in  several  cases 
demonstrated  poor  agreement  between  predicted  and  experimental  results.  Karlak 
performed  fracture  tests  on  a  series  of  quasi-isotropic  laminates  with  a  hole  of  different 
diameters  and  showed  that  the  characteristic  dimension  was  proportional  to  the  square  root 
of  the  hole  diameter.  Whitney  and  Kim  found  that  stacking  sequence  has  a  significant 
effect  on  the  notch  sensitivity  of  composite  laminates.  Pipes,  et  al  proposed  a  fracture 
model  based  on  the  principle  of  superposition  of  notched  strength  data  for  several 
composite  materials  and  laminate  stacking  sequences.  Through  the  development  of  a 
notched  strength-radius  superposition  method  and  the  definition  of  two  radius  shift 
parameters,  a  single  master  curve  was  developed  which  allows  the  prediction  of  the 
influence  of  laminate  stacking  sequence  upon  the  notched  strength  of  different  composite 
materials.  However,  this  approach  requires  additional  tests  to  determine  the  additional 
parameters  and  as  a  result  the  process  is  more  complicated  and  expensive.  The  accuracy  of 
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prediction  is  only  marginally  improved  for  those  cases  where  the  characteristic  dimension 
increases  with  notch  size. 

A  comprehensive  study  was  conducted  with  the  objective  of  collecting  a  massive 
set  of  experimental  results  and  comparing  the  data  with  predicted  results  generated  by 
several  commonly  used  fracture  models.  A  total  exceeding  2800  notched  strength  data  was 
collected  and  filed  in  a  computer  data  bank  for  graphite/epoxy,  boron/aluminium  and 
graphite  polyimide  laminates.  The  study  pointed  out  that  the  characteristic  dimensions  for 
different  material  systems  are  in  fact  independent  of  notch  size  in  most  cases  or  exhibit  only 
a  slight  increase  with  notch  size.  The  study  also  indicated  that  the  characteristic  dimension 
is  strongly  influenced  by  material  lay-up  patterns.  For  example,  the  values  of  the 
characteristic  dimension  of  ASl/3501-6  vary  from  199.6  mm  for  the  [Oj^,  lay-up  to  -0.06  mm 
for  the  [90]4  lay-up.  The  negative  value  of  the  characteristic  dimension,  which  has  no 
physical  meaning,  is  caused  by  the  condition  thata°°^/a„  <  1/K°^  ,.  This  condition  is  below 
the  lower  limit  of  Equation  4  and  therefore,  in  theory,  this  fracture  model  cannot  be  used 
to  predict  the  notched  strength  of  laminates  with  the  [90)4  lay-up.  However,  as  a  result  of 
the  curve-fitting  procedure  employed,  a  negative  value  was  obtained  which  gave  the  best 
correlation  with  the  test  data. 

Reviewing  the  summary  of  the  results  given  in  Reference  46,  it  is  observed  that  the 
effect  of  other  intrinsic  material  variables  on  the  notch  sensitivity  is  manifested  by  the 
significant  differences  among  the  values  of  the  characteristic  dimensions  for  identical  lay-ups 
obtained  from  different  sources,  suggesting  significant  influence  of  batch  variation  of 
constituent  properties,  fiber  volume  fraction,  fabrication  procedures,  etc.  For  example, 
Lagace*'^^'  obtained  a  value  of  1.4  mm  for  the  characteristic  dimension  of  ASl/3501-6  with 
a  lay-up  of  [0/±45]j.  which  differs  significantly  from  the  value  of  3.37  mm  obtained  by 
Cruse^'^**'  under  identical  test  conditions. 

In  this  work  the  effects  of  several  intrinsic  factors  were  controlled  by  adopting  a 
consistent  autoclave  fabrication  technique  and  determining  accurately  the  fiber  volume 
content.  The  quality  of  the  laminates  was  assured  by  using  ultrasonic  C-scan  inspection.  The 
inherent  scatter  of  the  test  data  was  accounted  for  by  using  adequate  sample  sizes  and  the 
Weibull  statistics  were  used  to  obtain  an  estimate  of  the  characteristic  dimension.  The 
established  fracture  model  was  assessed  by  using  it  to  predict  the  notched  strength  of  (a) 
IM6/5245C  laminates  containing  a  hole,  two  different  diameters  were  tested;  (b) 
IM6/5245C  laminates  containing  a  countersunk  hole:  (c)  laminates  fabricated  with  different 
second  generation  high  toughness  composites;  and  (d)  laminates  fabricated  with  a  first 
generation  composite,  '1300/5208. 

4.1  Notched  Strength  Prediction 

A  comparison  between  the  predicted  values  ofov^’"  /o„  and  the  experimental  results 
generated  for  IM6/5245C  laminates  is  shown  in  l  ig.  10.  I'he  predicted  \  allies  of  /a  , 
were  generated  using  Equation  4  for  a  range  of  hole  sizes.  'I  he  a^,  3.43  mm  established 
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earlier  in  Section  3  was  used  for  the  prediction.  This  a,,  was  obtained  by  corroborating  witli 
test  results  for  quasi-isotropic  laminates  containing  a  single  circular  hole.  Two  hole 
diameters,  6.35  mm  and  9.53  mm,  were  tested  for  this  purpose.  It  is  noted  that  the  predicted 
values  agree,  as  expected,  very  well  with  the  experimental  results  for  these  two  hole  sizes. 
However,  the  values  ofa^^/a^,  for  the  hole  diameters,  3.18  mm  and  12.7  mm,  are  predicted 
unconservatively.  The  discrepancy  between  prediction  and  experimental  results  appears  to 
be  related  to  the  small  scatter  of  test  data  resulting  from  small  sample  sizes  for  these  two 
hole  diameters.  The  number  of  tests  performed  for  each  of  these  two  diameters  is  three 
which  is  much  smaller  than  the  eleven  tests  performed  for  the  two  intermediate  hole 
diameters.  Although  confirmation  is  not  possible  without  additional  test  results,  it  appears 
that  the  small  scatter  may  not  be  representative  of  the  true  population  scatter  which  leads 
to  a  poor  correlation  with  the  experimental  results. 

The  fracture  model  with  a^^  =  3.43  mm  predicts  the  values  of  for  the 

countersunk  hole  very  well  when  the  countersunk  hole  is  approximated  by  a  uniform 
through-the-thickness  hole  with  a  diameter  equal  to  the  outer  diameter  of  the  countersunk 
hole.  This  finding  agrees  with  that  obtained  by  Nuismer  and  L;ibor 

The  results  for  predicting  the  values  of  o  /a  for  the  tw  o  '  -’cond  generation  high 
toughness  composites,  IM6/1806  and  1M6/F584.  and  for  a  first  generation  brittle  composite, 
T300/5208,  are  summarized  in  Fig.  11.  The  a„  =  3.43  mm  established  from  the  test  results 
of  IM6/5245C  is  used  in  the  fracture  model  for  predicting  the  values  of  for 

IM6/1806  and  IM6/F584.  The  a„  =  1.48  mm  established  from  the  test  results  of  AS4/3501-6 
is  used  for  predicting  the  values  of  a  /a  for  T300/5208.  As  shown  in  Fig.  11,  prediction 
agrees  with  experimental  results  only  in  the  case  for  1M6/F584,  a  composite  material  with 
mechanical  properties  similar  to  those  of  IM6/5245C  (Table  1). 

Using  the  procedure  discribed  in  Section  3.  the  characteristic  dimensions  for 
IM6/1806.  IM6/F584  and  3300/5208  were  determined  to  be  4.9  mm.  3.3  mm  and  2.7  mm 
respectic Icy.  It  must  be  emphasized  that  these  values  ha\e  been  determined  based  on  3  to 
4  tests  for  each  material  and  they  may  not  be  strong  estimates  of  the  true  a  for  these 
materials.  Using  these  values  the  prediction  of  was  performed  over  a  range  of  hole 

sizes  for  these  materials  and  the  results  arc  presented  in  Fig.  12.  It  is  observed  that  a 
decrease  in  a„  results  in  an  increase  in  notch  sensitivity,  'riicrcfore  a  ranking  of  the  notch 
sensitivity  of  composite  materials  based  on  a„  is  possible. 

Although  the  a,,  of  a  composite  material  can  be  correlated  to  notch  sensiticity,  it 
cannot  indicate  the  change  in  failure  mode  resulting  from  a  change  m  notch  sensitivity.  .\s 
shown  in  big.  13  an  increase  in  notch  sensitivity  resulting  from  an  increase  in  hole  diameter 
from  3.18  mm  to  12.7  mm  corresponds  to  a  change  in  failure  mode  from  slant  fracture  to 
flat  fracture.  l  or  intermediate  hole  sizc>,  mixed  modes  involving  slant  and  flat  fractures 
predominate.  As  shown  previously  a  constant  value  of  a„  =  3.43  mm  can  be  used  to  predict, 
with  mixed  results,  the  notched  strength  for  all  these  hole  sizes.  The  countersinking  of  a  hole 
with  a  nominal  diameter  =  6.71  mm  results  in  a  significant  increase  in  notch  ,^ensitivitv  and 
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a  flat  fracture  mode  as  shown  in  Fig.  14.  The  same  value  of  a^,  =  3.43  mm  can  be  used  to 
obtain  a  good  correlation  with  the  lest  results  by  approximating  the  countersunk  hole  by  a 
uniform  through-the-thickness  hole  with  a  diameter  equal  to  the  outer  countersink  diameter. 


5.0  CONCLUSIONS 


An  investigation  was  conducted  which  led  to  the  determination  of  the  two  parameters 
of  a  fracture  model  for  the  prediction  of  the  tensile  notched  strength  of  several  graphite 
fibre  reinforced  composites.  Key  conclusions  drawn  from  this  work  are: 

1.  The  merit  of  the  two  parameter  fracture  model  based  on  the  average  stress 
criterion  is  the  inherent  simplicity  of  the  model.  It  was  found  to  be  capable 
of  predicting  the  tensile  strength  of  composite  laminates  containing  notches 
of  various  sizes. 

2.  The  two  parameters,  which  are  the  unnotched  strength  of  the  subject  laminate 
and  the  characteristic  dimension  of  the  material,  must  be  determined  before 
the  model  can  be  used  for  prediction. 

3.  The  determination  of  the  characteristic  dimension  requires  a  data  base 
including  the  notched  strength  and  elastic  ply  properties  of  the  subject 
laminate.  Since  the  characteristic  dimension  is  dependent  on  a  variety  of 
intrinsic  and  extrinsic  factors,  it  must  be  determined  by  experiment  for  each 
new  composite  material,  laminate  configuration,  etc. 

4.  A  statistical  procedure  based  on  the  two  parameter  Wcibull  distribution  was 
used  to  obtain  an  estimate  of  the  characteristic  strength  parameters  for  the 
notched  and  unnotched  strengths  of  the  subject  laminate.  The  ratio  of  the 
characteristic  strength  parameters  was  used  to  calculate  the  characteristic 
dimension. 

5.  The  characteristic  dimension  can  be  used  as  an  indicator  for  the  notch 
sensitivity  of  different  composite  materials.  A  decrease  in  the  characteristic 
dimension  was  found  to  indicate  an  increase  in  notch  sensitivity. 

6.  The  second  generation  high  toughness  composites  used  in  this  study  which 
included  the  IM6/5245C,  1M6/F584  and  IM6/1806  were  significantly  less 
notch  sensitive  than  the  first  generation  composite  T300/5208.  The 
characteristic  dimensions  of  these  composite  materials  were  larger  than  that 
of  the  T300/5208. 
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7.  Approximating  the  countersunk  hole  by  a  uniform  through-the-thickness  hole 
with  a  diameter  equal  to  the  outer  diameter  of  the  countersunk  hole  allows 
the  fracture  model  to  be  used  for  predicting  the  tensile  notched  strength  of 
a  countersunk  hole. 
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TABLE  1.  MECHANICAL  PROPERTIES  OF  COMPOSITE  MATERIALS 


(a)  Elastic  Moduli 


Material 

Name 

0°  Tensile 
Modulus 
(ETi),GPa 

90°  Tensile 
Modulus 
(Et2),  GPa 

Shear 

Modulus 

(G,2),GPa 

Poisson 

Ratio 

(M-12) 

Narmco 

IM6/5245C 

166.2 

8.3 

5.5 

0.31 

Hexcel 

IM6/F584 

171.0 

9.0 

5.9 

0.32 

Cyanamid 

IM6/1806 

155.1 

10.3 

3.9 

0.32 

Hercules 

AS4/3501-6 

140.0 

8.2 

6.2 

0.30 

(b)  Ultimate  Strengths 


Material 

Name 

0°  Tensile 
Strength 
(Xti),  MPa 

0°  Compressive 
Strength 
(Xci),MPa 

90’’"  Tensile 
Strength 
{Xt2),  MPa 

90^  Compressive 
Strength 
(Xc2),  MPa 

Inplane  Shear 
Strength 
(X12K  MPa 

Narmco 

1M6/5245C 

2,620.0 

1,278.3 

60.0 

220.6 

1 17.9 

Hexcel 

1M6/F584 

2,546.2 

1.340.0 

47.6 

232.4 

122.7 

Cyanamid 

IM6/1806 

1,840.9 

1,178.3 

40.7 

199.3 

91.0 

Hercules 

AS4/3501-6 

2,144.3 

824.6 

46.2 

172.4 

1 10.3 

7 


TABLE  2.  STRESS  CONCENTRATION  FACTORS  FOR  COMPOSITE  LAMINATES 
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Lay-up  II  =  [+45/0/-45/()/90/0/+45/0/-45/0]  2s 


TABLE  3.  UNNOTCHED  TENSILE  TEST  RESULTS 


Material 

Name 

Lay-Up* 

Specimen 

Number 

Ultimate  Tensile 
Strength 
(Co),  MPa 

Modulus 
@  .004  Strain 
(Ey),  GPa 

IM6/5245C 

I 

2 12- A 

846.7 

N/A 

IM6/5245C 

I 

212-B 

853.6 

N/A 

IM6/5245C 

I 

212-C 

842.5 

N/A 

IM6/5245C 

I 

212-D 

766.0 

N/A 

IM6/5245C 

I 

212-E 

786.7 

N/A 

IM6/5245C 

I 

212-F 

862.5 

N/A 

IM6/5245C 

I 

212-G 

787.4 

N/A 

IM6/5245C 

I 

212-H 

787.4 

57.9 

IM6/5245C 

I 

211-X 

836.3 

N/A 

IM6/5245C 

I 

211-Y 

912.2 

IM6/5245C 

I 

211-Z 

815.6 

IM6/5245C 

11 

1,274.2 

90.3 

IM6/5245C 

II 

1,370.7 

89.6 

IM6/5245C 

II 

157-Y 

1,200.4 

91.0 

IM6/5245C 

II 

157-Z 

1,256.2 

88.3 

AS4/3501-6 

II 

233-A 

1,150.0 

78.6 

AS4/3501-6 

II 

233-B 

1,179.0 

83.4 

AS4/3501-6 

II 

233-C 

1,079.0 

83.4 

AS4/3501-6 

II 

233-D 

1,170.0 

80.7 

AS4/3501-6 

II 

233 -E 

1,116.3 

78.6 

AS4/3501-6 

II 

233-F 

1,137.6 

80.7 

AS4/3501-6 

II 

233-G 

1,155.6 

81.4 

AS4/3501-6 

II 

233-H 

1,261.1 

82.0 

AS4/3501-6 

II 

233-1 

1,150.0 

79.3 

AS4/3501-6 

II 

233-J 

1,163.8 

81.4 

AS4/3501-6 

II 

233-K 

1,119.7 

79.3 

AS4/3501-6 

II 

233-L 

1,230.0 

N/A 

*  see  Table  2 
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TABLE  4.  TEST  MATRIX  FOR  NOTCHED  SPECIMENS 


Material 

Name 

Lay-Up* 

Hole  Type 

Diameter 

mm 

Number 
of  Tests 

IM6/5245C 

I 

Circular 

3.18 

3 

IM6/5245C 

I 

Circular 

6.35 

11 

IM6/5245C 

I 

Circular 

9.53 

11 

IM6/5245C 

I 

Circular 

12.7 

3 

IM6/5245C 

II 

Countersunk 

100° 

6.71  (min.) 
13.46  (max.) 

9 

IM6/1806 

I 

Circular 

6.35 

3 

IM6/F584 

I 

Circular 

6.35 

4 

T300/5208 

I 

Circular 

6.35 

3 

AS4/3501-6 

II 

Circular 

6.71 

11 

*  see  Table  2 
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TABLE  5.  NOTCHED  TENSILE  TEST  RESULTS 
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TABLE  6.  RESULTS  OF  THE  DETERMINATION  OF  THE  WEIBULL  PARAMETERS  FOR  UNNOTCHED  TEST  DATA 
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TABLE  7.  RESULTS  OF  THE  DETERMINATION  OF  THE  WEIBULL  PARAMETERS  FOR  NOTCHED  TEST  DATA 

(a)  MATERIAL  =  IM6/5245C 
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see  Table  2 
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TABLE  8.  CHARACTERISTIC  DIMENSIONS  FOR  IM6/5245C  AND  AS4/3501-6 


Material 

Name 

Lay-Up* 

Diameter 

mm 

Characteristic 
Strength  Ratio 

fc/pj 

Characteristic 
Dimension 
(Ao),  mm 

IM6/5245C 

I 

6.35 

0.60 

3.43 

IM6/5245C 

I 

9.53 

0.54 

3.43 

AS4/3501-6 

II 

6.71 

0.46 

1.48 

*  see  Table  2 
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FIG.  1 :  AVERAGE  STRESS  FAILURE  CRITERION 
FOR  A  CIRCULAR  THROUGH-THE-THICKNESS  HOLE 
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STRENGTH  RATIO 


FIG.  2:  EFFECT  OF  THE  CHARACTERISTIC  DIMENSION,  a^,  ON 
NOTCHED  STRENGTH  PREDICTION  FOR  LAMINATES  CONTAINING 

oo  oo 

A  CIRCULAR  HOLE  HAVING:  (a)  Kj  =  3.0;  (b)  =  6.0 


SPECIMEN  #  233A 
AS4/3501-6 

[+45/0/-45/0/90/0/-t-4S/0/-45/0] 
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FIG.  3:  A  TYPICAL  STRESS  vs.  STRAIN  PLOT  FOR  AN  UNNOTCHED  TENSILE  TEST 


FIG.  4:  A  TYPICAL  SETUP  FOR  A  NOTCHED  TENSILE  TEST 


[+45/0/-45/90]  6S 
►  TENSILE  TEST  DATA 


WEIBULL  PROBABILITY  PLOT  FOR  IM6/5245  UNNOTCHED  STRENGTH  DATA 
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FIG.  8:  WEIBULL  PROBABILITY  PLOT  FOR  IM6/5245C  NOTCHED  STRENGTH  DATA 


AS4/3501-6 

[+45/0/-45/0/90/0/+45/0/-45/0]  2S 
¥  TENSILE  TEST  DATA 


I 


(BdN  NT 

oo 
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FIG.  9:  WEIBULL  PROBABILITY  PLOT  FOR  AS4/3501-6  NOTCHED  STRENGTH  DATA 


DATA  RANGE 
CIRCULAR  HOLE 


FIG.  10:  COMPARISON  BETWEEN  PREDICTED  ^  AND  TEST  RESULTS  FOR  IM6/5245C 
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FIG.  12:  EFFECTS  OF  CHARACTERISTIC  DIMENSION  ON  NOTCH  SENSITIVITY 


FIG.  13:  EFFECTS  ON  FAILURE  MODE  RESULTING 
FROM  INCREASING  HOLE  SIZE 
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FIG.  14:  FLAT  FRACTURE  OF  AN  IM6/5245C  LAMINATE 
CONTAINING  A  COUNTERSUNK  HOLE 
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SUMMARY/SOMMAIRE 


This  report  describes  the  determination  of  empirical  parameters  for  a  fracture  model  which  can 
be  used  to  predict  the  uniaxial  tensile  stren^h  of  composite  laminates  containing  notches  of 
various  sizes.  The  parameters  required  are  the  characteristic  dimension  and  the  unnotched 
tensile  strength  of  the  subject  laminate.  A  data  base  including  notched  and  unnotched  strengths 
of  the  subject  laminates  and  elastic  ply  properties  was  generated  for  the  determination  of  these 
parameters  using  a  statistical  procedure  based  on  the  two  parameter  Weibull  distribution.  Three 
second  generation  high  toughness  graphite  fibre  reinforced  composites,  1M6/5245C,  IM6/F584 
AND  IM6/1806,  and  two  first  generation  brittle  graphite/epoxy  composites,  AS4/3501-6  and 
T300/5208,  were  used  to  fabricate  the  specimens  containing  either  a  through-the-thickness 
cylindrical  hole  or  a  countersunk  hole. 
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